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We extend and apply a previously developed method for a semiclassical treatment of a system with 
large spin 5*. A multisite Heisenberg Hamiltonian is transformed into an effective classical Hamilton 
function which can be treated by standard methods for classical systems. Quantum effects enter 
in form of multispin interactions in the Hamilton function. The latter is written in the form of 
an expansion in powers of J/{TS), where J is the coupling constant. Main ingredients of our 
method are spin coherent states and cumulants. Rules and diagrams are derived for computing 
cumulants of groups of operators entering the Hamiltonian. The theory is illustrated by calculating 
the quantum corrections to the free energy of a Heisenberg chain which were previously computed 
by a Wigner-Kirkwood expansion. 

PACS numbers: 05.30.-d, 75.10.Hk, 75.10.Jm 



I. INTRODUCTION 

Quantum effects in systems with localized spins de- 
pend on the size of the spin S. They are strongest for 
5 = 1/2 while systems behave classically in the limit 
S ^ CO. There are, however, numerous compounds hav- 
ing intermediate or high spin values, which we would like 
to treat semiclassically, i.e., classically with quantum cor- 
rections expanded in powers of 1/5. Different proposals 
have been made in the past how this can be achieved 
We have recently suggested a method by which 
we can derive effective classical Hamilton functions for 
spin Hamiltonians by using spin coherent states and cu- 
mulants 1^. The quantum partition function is thereby 
replaced by a classical one, and the trace over a complete 
set of quantum states is replaced by an integration ober 
the classical spin vectors. The mathematical structure 
of these two objects is quite different, and the transition 
from one to the other was investigated by Lieb [Q who 
obtained upper and lower bounds for the quantum par- 
tition fimction for large 5. He did not consider the 1/5 
corrections to the classical partition function which arc 
discussed here. 

The previous paper |^ dealt with one- and two-spin 
systems and their semiclassical description when 5 is 
large. Here we extend these calculations to a many-spin 
system. The necessary formalis m is described in Sec. 
y. The method is applied in Sec. Ill to the Heisenberg 
Hamiltonian in an applied magnetic field first stating the 
rules for evaluating the required cumulants. As we will 
show this can be done by means of a diagrammatic tech- 
nique. The explicit form of the corresponding classical 



Hamilton function is derived in Sec. IV. All terms of 



order J[J/{TS)]^/S" with n = 0, 1, 2 are included in the 
effective Hamilton function. Obviously, the semiclassical 
description of large-spin systems breaks down in the low- 



temperature limit, where quantization of energy levels 
becomes important. To check the method, we use in Sec. 
VI our effective Hamilton function in order to calculate 
the free energy and other thermodynamic properties of a 
spin chain to order 1/5^. These results were previously 
computed by a Wigner-Kirkwood expansi on w ithout us- 
ing effective classical Hamiltonians. Sec. VII contains a 
discussion and the conclusions. 



II. SPIN COHERENT STATES AND 
CUMULANTS 

To achieve our goal, we use two theoretical tools. The 
first one is the coherent-state representation of quantum- 
statistical averages. On each lattice site i, we introduce 
spin coherent states jrii), i.e., states with the maximal 
projection on the axis pointing in the direction of the 
unit vector n^. The direct product of these states |{ni}) 
approaches the "classical" state of the spin system in the 
limit 5 oo. On the other hand, the basis of coherent 
states is complete (and even overcomplete), so that it 
contains all quantums states. The quantum-statistical 
averages of the system can be rewritten in the coherent- 
state basis with the help of the unity operator 



25+1 
An 



dn |n)(n| 



(1) 



For a single-spin system, the trace of an operator A over 
any complete orthonormal basis \m) becomes 



trA^ '^{m\A\m) = 

m 

25 + 1 



25+1 
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dn ^^(m|A|n) (n|TO) 



(2) 
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Therefore the partition function for a many-spin quan- 
tum Hamiltonian H can be written as 

^ ^ \r^) I n^"' ({n,}|exp(-/3iJ)|{n,}), (3) 

where (3 = 1/T. It has the same form as the one for clas- 
sical systems, provided one defines the effective classical 
Hamilton function Ti by the relation 

m - -ln({nj|exp(-/3i?)|{n,}} 

= ({n,}|l-exp(-/3i7)|{njr (4) 

Here the superscript c denotes the cumulant of a matrix 
element (see below). Cumulants constitute the second 
theoretical tool used in this paper. 

The above formula is exact for any quantum spin sys- 
tem. In the limit S* — > oo, reduces to the usual clas- 
sical Hamiltonian. A mathematically rigorous treatment 
of this limiting transition can be found in Lieb's paper, 
Ref. Q Here we will obtain quantum corrections to the 
classical Hamiltonian using a systematic 1/S expansion 
for S ^ 1. In the same context, a technique using partial 
differential equations for the density matrix was applied 
in Ref. [Q. But an explicit derivative- free expression for 
the effective classical Hamiltonian was obtained in Ref. 

only for the spin-chain model under the condition that 
neigbouring spins are nearly coUinear. The method pre- 
sented here does not use derivatives and yields explicit 
results for any lattice. 

The effective Hamilton function is obtained in form of 
a cumulant expansion by expanding in Eq. (Q) the op- 
erator exp(— /3i?). It is worth noting a basic relation 
between semiclassical theory. A characteristic property 
of quantum mechanics is that the expectation value of a 
product AiA2 of two observables with respect to a quan- 
tum state is generally distinct from the product of 
the expectation values of Ai and A2. The difference is 
just the cumulant, i.e., 

(AiA2)^-(AiA2)-(Ai)(A2)^0. (5) 

In the classical limit this difference vanishes. Therefore 
we expect that a theory formulated in terms of cumulants 
is particularly suitable for a semiclassical expansion. 

Let us recall the definition and basic properties of cu- 
mulants. Averages or matrix elements (. . .) can be ex- 
pressed through cumulants {. . -Y as follows 

{A) = {AY, = {A^A^Y + {Ai){A2), 

{A^A^A^) = {AiA2A^Y + {Ai){A2A:iY + {A2){A^A^Y 
+ {A^){A^A2Y + {A,){A2){A^), (6) 

etc., where Ai are classical stochastic variables or 
quantum-mechanical operators. The averaging above 
is performed over a classical distribution function or 
weighed over quantum states. A more detailed discus- 
sion of cumulants can be found in refs. ||^-^. Cumulants 
can be obtained by differentiation of a generating func- 
tion, i.e., from 



at Ai = . . . = A„ = 0, in contrast to averages (...), which 
are given by a similar expression without logarithm. By 
multiplying by 11"=! '^'^'/^i- ^^'^ summing over all ki = 
0, . . . , 00 without the term fci = . . . = fc„ = one obtains 
the identity 

(e^i^i . . . e^"^" - 1)'= = ln(e^i^i . . . e^"^"). (8) 

The second line of Eq. is a particular case of this 
formula. 

Let us consider cumulants of spin operators with re- 
spect to spin coherent states. These cumulants have 
an especially simple form when the spin-operator com- 
ponents Sz (axis z along n) and S± = Sx i iSy are 
used. A cumulant vanishes if it is of the form (. . . 5'+)'^, 
(. . . S+S-)'', (. . . Sz)'', etc., where . . . stands for any com- 
bination of spin operators. That is, the number of S+ and 
S- in the cumulant must be balanced in order to give 
a nonzero result. If, however, this balance is achieved 
already within a subgroup of operators on the right or 
left side of the operator list, the cumulant vanishes, too. 
Non-vanishing cumulants are 

{Sz) - s, {s+s:s-Y = 25(-l)", 

{S+S^S+Sl^S^S'/S-Y = -4S'(-l)"+"'(-2)'", 

{s+s+s+S-S-S-Y = i{s+s+S-S+S-S-Y = 245, (9) 

etc. The above results can be obtained recurrently, using 
Eq. (p|) and the commutation relations [S'^,.?-] = —S- 
and p+,S'_] = 2Sz, which remain valid inside cumu- 
lants. For scaled spins S/S, each nonvanishing cumulant 
containing n spin operators scales like l/S""^^. 

An expansion of Ti. in Eq. (^) in powers of 1/S is ob- 
tained from the Taylor scries 

I3H = {1- cxp{-PH)Y 

= mY - ^{HHY + ^{HHHY + ■■■, (10) 

where we have introduced the shorthand notation (...) = 
({rii}! . . . |{ni}). Here the first term on the right-hand 
side is the classical energy of the spin. As seen from Eq. 
(^, increasing powers of IJS appear in each order of the 
expansion. In fact, Eq. (Ilfl) is an expansion in powers 
of J/{TS) and it breaks down at low temperatures. For 
5 S> 1, the range of convergence is much larger than 
that of the high-temperature series expansion. For the 
calculation of thermodynamic quantities one can further 
expand exp(— /37i) in powers of 1/ S. In that case the sta- 
tistical weights of different spin orientations are described 
by a purely classical Hamiltonian, whereas quantum ef- 
fects manifest themselves in corrections to the density of 
states. 
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III. THE CUMULANT EXPANSION 

Let us consider a spin Hamiltonian of the Heisenberg 
form 



H — ^ Hi - Si — - ^ Jij Si • Sj . 



(11) 



To implement the cumulant expansion of Eq. (|1C|), it is 
convenient to express the spin operator on each site i in 
the coordinate system with the z axis along the coherent- 
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FIG. 1. Diagrammatic representation of the cumulant 
expansion for the quantum Heisenberg magnet. Wavy lines 
stand for the exchange interaction Jij, small solid circles are 
cumulants of spin operators. 



Sj = ^ riia-SiCi n± = (ria; T i%)/2, (12) 

where and ar e ap propriate transverse basis vectors. 
Insertion into Eq. ( J10| ) leads to expressions of the type 



2! 22 



ii'jj' OiiOi^l OLjOL^I 

X {{SiaiSi'a^,){Sjaj Sj'aji)) 7 (13) 



as illustrated by the second-order pure-exchange term. 
The brackets inside the cumulant imply that this cu- 
mulant is defined with respect to the two pairs of the 
spin-operator components and not with respect to four 
single operators. To calculate such cumulants of com- 
posite operators, it is convenient to (i) express them 
through ordinary matrix elements [see, e.g., Eq. (j^)] and 
(ii) express matrix elements through cumulants of single 
spin-operator components [see Eqs. (^]. The latter are 
nonzero only if all the operators in the cumulant belong 
to the same site in which case they are readily given by 
Eqs. (^). One can see that in several lowest orders of 
the cumulant expansion, the summation over the spin- 
component indices a in Eq. ( p^ reduces to a single real- 
ization which gives a nonzero result. The result of such 
a procedure applied to Eq. (fsl) is the following: 



22 



+ ^{AA,r{A,A,,r, 



(14) 



where Ai = Sia^ . Hereby we have taken into account 
that operators Ai and Aii belong to different sites so that 
terms of the type {AiAi')'' vanish. On the rhs of this ex- 
pression, one could add terms differing by permutations 
of the indices and/or jjf- These terms make the 
same contributions to Eq. (^^ as those present in Eq. 
(p^); Instead of writin g th em explicitly, they lead to the 
prefactor 1/22 (|lj)- Generally, the possibility of 

permutating spin operators connected by the exchange 
interaction effectively cancels the coefficient 1/2 in Eq. 
(|n|). The factor 1/2 in front of the last term of Eq. (|l|) 
appears because of the permutation of both i, i' and j, j' 
which does not generate a new term of this type. For the 
cumulant with three groups of spin operators one obtains 



1 

23 



{{A,A,){A,A,,)iAiAi,)Y 



{A,){A,A,nA,,A,y{A,) 
+ {Ai){A'Air{A,Ai,r{A,,) 
+ {A,){AA,'y{A,Air{A,) 
+ {AAjAir{A,){Aj,){Ai,) 

+ {AA,y{A,AyAiy{Ai,) 

+ {A,Aiy{A,A,Ai,y{Aj,) 
+{AjAir{AarAi,r{A,) 
+ {A-AjriAyAiyiAAur 

+ ^{A,A,Air{AeA,,A,r. 



(15) 



The rhs of Eqs. (14) and ( |15| ) are constructed according 
to a principle which can be formulated in a diagrammatic 
language: The operators A are either contracted into 
cumulants or connected by the interaction lines. Sim- 
ilarly to other diagram techniques, there are no terms 
consisting of disconnected parts, e.g., there is no term 
{A,){A,,){Aj){Aj,) in Eq. (]l|). Note that the order of 
the operators in the cumulants on the rhs of Eqs. ( p^ ) 
and ([l5| ) is the same as on the left-hand side (Ihs). This 
fact enables us to write down immediately equations of 
the form of Eqs. ( [T^ ) or (^5|), without explicitly perform- 
ing steps (i) and (ii). The diagrammatic representation 
of the cumulant expansion of Eq. ( p^ ) for the quantum 
Heisenberg magnet in the zero-field case is shown in Fig. 
|l|. Diagram is the first order of the cumulant expan- 
sion in H. It will be shown below that this diagram yields 
the zeroth order inl/S for the Hamilton function Ti., i.e., 
the classical Hamiltion function Ti^^K Diagrams la and 
lb correspond to the two terms in Eq. (M). Diagram 2a 
represent the three first terms in Eq. ([l5|). These terms 
have the same topology and differ only by the order of 
spin operators. This reflects the quantum nature of the 
latter. Diagram 2b corresponds to the fourth term in Eq. 
(|l5|). Diagram 2c represent the fifth, sixth, and seventh 
terms in Eq. ( p^ ) . Diagrams 2c? and 2e represent the two 
last terms, respectively. 

Cumulants corresponding to the mixed field-exchange 
terms in the cumulant expansion can be treated in a sim- 
ilar manner. One obtains 



1 



{{AMA,r = {A,){A,A,Y 
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J.- uJ> LJ 



JH 



JHH-1 



.1 ' 

JHH-2 JJH-4 



JJH-1 



JJH-2 




applying some vector algebra (for details see Appendix). 
This is not surprising since n is the only vector specifying 
the spin coherent state. Thus one can rewrite in the 
form 



AS 



FIG. 2. The same as in Fig. ^ for the mixed field-exchange 
terms. Straight lines represent an external magnetic field. 



-{A,{A,A,)Y ^ {A,){A,A,Y 

]^{{A,A,,)AjAiY = {A,){A,,A,AiY + {A,AiY{A,,A,Y 

^{{MA,){A,A,.)AiY = {AAiY{A,A,Y{Ar) 

+ {A,){A,A,AiY{A,,) + {A,){A,A,Y{AyAiY 

+ {A,A,AiY{A,AyY, (16) 

etc. The diagrams corresponding to the rhs of these ex- 
pressions are shown in Fig. ^j. For the last equation of 
Eqs. (p^), in particular, the diagram corresponding to the 
nth term {n — 1, 2, 3, 4) on the rhs is labeled by JJH-n. 



IV. THE EFFECTIVE HAMILTONIAN 

Let us now proceed to the calculation of the different 
terms in the cumulant expansion of Eq. ([l0| ) . The term 
of the first order in H yields the classical Hamiltonian 

7^(0) = _^h,.n,-i^J,,n,.n„ (17) 



which is the zeroth order of the 1/S expansion. Here 

hi = 5Hi, Jij = S^Jij (18) 

are the reduced magnetic field and the exchange interac- 
tion, respectively. In second and higher orders in H of 
the cumulant expansion, 1/5' corrections appear. With 
the help of Eqs. (|l|), (0), and (§) one obtains 

/3 



/3 



- Jy [(h, • ni+)(n,_ • Uj) + (h, • n,_)(ni+ • n^)] 

- ^ J^JJJl{n, ■ nj+)(nj_ • n;) 

ijl 

52 



(19) 



In this expression, as well as in all other expressions of 
this type, the transverse components of the coherent- 
state vectors, n±, can be expressed in terms of n after 



45 



45 



J2 ^Kh* • nj) - (hj • ni)(n, • n^) 

ij 

+ (hj • rii) - (hj ■nj){nj ■ n,)] 



ijl 



f3 



1652 



5]j2[l-(n,-n,)]2 



(20) 



One can see that the effective classical Hamilton func- 
tion corresponding to the Heisenberg quantum Hamilto- 
nian, Eq. (pT|), is of a non-Heisenberg form. In particular, 
many-spin interactions appear. 

Let us now proceed to the third order terms of the 
cumulant expansion. At first we consider the zero-field 
case. Equation (|l3) generates terms of different topology 
in 7^(2) (geg Fig^) 

Ti.^^^ = -— ^ Jy Jj7Jir!$2a(ni,nj,n;,n„) 



ijln 



f3' 



1252 



Ji7iJjnJln^2b{T^i,'nj,ni,nn) 



ijln 



2453 



4853 



4854 



^ Jy Jj/$2c(ni,nj,ni) 



ijl 



Y JijJjiJii^2d{ni,nj,ni) 



ijl 



^ Jf^$2e(n»,nj). 



(21) 



Here the 2a term is due to the first, second, and third 
terms of Eq. (15), the 2b term is due to the fourth term 
of Eq. (|l^), the 2c term is due to the fifth, sixth, and 
seventh terms of Eq. (|l5|), the 2d term is due to the eighth 
term, and the 2e term is due to the nineth term. For the 
functions $ in Eq. (|2l|) we find 



16 

y 



(rii • nj+)(nj_ ■ n;+)(n/_ • n„) 
+ (rii • nj+)(nj_ • n/_)(n;+ • n„) 
-I- (rij • nj_)(nj+ • n;+)(n(_ • n„) 



: (rii • n„) - (n,; • n;)(n; • n„) 
(rii • nj){nj ■ n„) + (n^ • nj){nj ■ n;)(n; • n„) 

^ [(n^ ■ ni){nj ■ n„) - (n^ ■ n„)(nj • n;)] (22) 
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= -4(nj ■ n„+)(nj • n„)(n; • n„_) 
= (n^ ■ n„)(nj ■ n„)(n, • n„) - 



(rij • n„)(nj • n;) 



{uj ■ n„)(ni • n;) + (n; • n„)(ni • n^) 



(23) 



^2h,a = (hi • n;) - (hj • nj){nj ■ n;) 
+ (hi • Ui) - {hi ■ nj){nj ■ n^) 

- [(hi • rij) + (hi • n;)] [(n^ • n;) - (n^ • nj){nj ■ n;)] 
(rii • n;)[(hi • rij) + (h/ • rij) 



$2c = -16[(ni+ ■ nj+)(ni_ • nj)(nj_ • n;) 

+ (ni_ • nj_)(ni+ • nj)(nj+ • n/) 

+ (ni+ • nj+)(ni_ • nj-)(nj • n/)] 
= [1 - (n, • n,)]{[3(n, • n,) - l](n, • n,) - 2(n, • n,)}, (24) 

$2d = 4^(ni+ • nj+)(nj_ • n;+)(n(_ • ni_) 

= (n, ■ Ujf + (Uj -11;)^ + (n; • Uif 

- (rij • nj){nj ■ ni){ni ■ n,;) 

(rii • rij) + (rij ■ n;) + (n; • n^) + (n.j ■ nj){nj ■ ni) 



{h, ■ ni){n, ■ Uj) - {hi ■ ni)(n; ■ 

2 " 



2(ni • nj){hj ■ nj){nj ■ n;) 



(tii •n,)(hj -rij) 



(rij • n;)(n; • n^) + (n; • ni)(ni • n^) 



(25) 



and 



$2e = 16(ni+ • nj+)(nj • nj)(ni_ ■ nj_) 
= (n,.n,)[l-(n..n,)]2. (26) 

The coefficient can be rewritten in the form 
$2d = - [niWjnif + [1 - (ni • nj)][l - (n^ ■ n;)(n; • rii)], 

(27) 

where [abc] = a • (b x c) is the mixed product and an 
appropriate symmetryzation is impHed. 

The field-dependent terms of order , JH^ and, J^H 
m can be calculated with the help of Eq. ( |T6| ) . The 
result has the form 



H 



12S'2 



^(rii • hj)[n,, X hj]2 



~ ^ H *2/i2 (n,, nj, h„ hj) 



^ X! '^!j'^ji*2/i,a(n^, rij, n;, h,, h^, h; 



ij; 



4853 



(28) 



where 



$2^2 = (n^ • rij) [(h^ • n^)^ + {hj ■ n^)^] 
2 

- - [(h, • n,)(h, -11^) + (hj • nj)(hj • n,)] 
-i(n,.n,)[/i? + 

+ (hi ■ hj) - (hj • ni)(ni ■ h^) - (hi • nj)(nj • hj) 
+ (hi ■ ni){n, ■ nj){nj ■ hj) 

+ i [(h, • nj)(hj • n,) - (h, • hj){n, ■ n^)] , (29) 



+ (nj • ni){hj ■ n,) + (rij • n,)(hj ■ n,) 



(30) 



and 



<^2h,b = -[1 - (rii • nj)]{2[(hi -11^) + (h^ • m)] 

+ [1 - 3(n, • n,)][(h, • n,) + (h, • n,)]}. (31) 

For a homogeneous field, hi = h, the last expression 
becomes 



2h.b 



= -3[h.(n, + n,)][l-(n,.n,)]' 



(32) 



As was said after Eq. (|l^), cumulant expansion is, in 
general, an expansion in P/S. In particular, in zero field 
in the nth order of the cumulant expansion, terms of or- 
der J{j3J /S)^~^ appear. On the other hand, there are 
also terms carrying additional powers of 1/ S, such as the 
last term of Eq. ( pO| ) and the three last terms of Eq. ([2]]). 
This feature results, formally, from reexpressing of cu- 
mulants of composite operators by ordinary cumulants, 
Eqs. ( [T^ ) and ([Tsl), which is specific to many-spin sys- 
tems. For 3> 1, the terms containing powers of (3J/S 
without additional factors 1/S dominate over the other 
terms in the whole range of temperatures. Those terms 
are given by "tree" diagrams, such as 0, la, 2a, and 2b in 
Fig. These diagrams are maximally branched and they 
do not contain loops or parallel interaction lines. Non- 
tree diagrams can be obtained from the tree diagrams by 
joining two solid circles into one. This leads each time to 
appearance of a spin cumulant of a higher order and thus 
to an additional factor of l/S" [see Eq. (^]. Minimally 
branched diagrams are those consisting of several paral- 
lel wavy lines, such as the diagrams lb and 2e in Fig. 
These diagrams result in terms of order J{I3J / S"^)"^^^ in 

n. 



V. PHYSICAL QUANTITIES 

Quantum-statistical averages of operators describing 
various physical quantities can be obtained by differen- 
tiation of the partition function Z or its logarithm with 
respect to appropriate parameters. For the internal en- 
ergy U = {{H)) one obtains 

U = -dlnZ/dl3= (n*), (33) 

where (...) denotes a classical thermal average and 
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(34) 



differs from H = H^^^ +H^^'> +H'^^^ + ... The scaled mag- 
netization per site, m = {{S))/{SN) = Ej((S»))/(5'A^), 
is given by 



1 91nZ 1 , 



n = (35) 



where h is the homogeneous part of the magnetic field 
defined formally by => h.; + h. Here, n* is not just 
n = rii but contains quantum corrections from all 
orders of the cumulant expansion. In particular, to first 
order of the cumulant expansion, one obtains from Eqs. 
(0) and ^) 



(36) 



The reduced correlation function of different spin com- 
ponents on different lattice sites can be written as 



dH on 



dhia dhjf3 dhiadhjfj 



(37) 



One notices that in order to calculate a correlation func- 
tion, it is insufficient to perform a classical thermal av- 
erage of niatijp or even of n*^n*p. The last term of Eq. 
( ^ ) makes a contribution to third order in the cumu- 
lant expansion due to the terms of the type JH^ [see Eq. 



VI. APPLICATION TO THE SPIN CHAIN 

The isotropic spin chain in zero magnetic field is a sim- 
ply solvable model in the classical limit |^ . The effective 
quasiclassical Hamiltonian discussed in this paper can be 
used to analytically calculate 1/ S corrections to the clas- 
sical results. The 1/ S expansion of the partition function 
has the form 



(38) 



where the averages are performed with respect to the 
classical Hamiltonian 



2o 



25" +1 
47r 



N 



^0, 



(39) 



and Zq is the partition function of the classical sys- 
tem. For the open spin chain the latter is given by 
Zo = 47r(47rsinh(^)/C)^-i with ^ = /3J. To order 
one should use for the linear W^^^ term in Eq. (pSl) the 



third and fourth terms of Eq. (po| ), for the quadratic TY^^^ 
term the third term of Eq. (^0|), and for the linear Ti'^^ 
term the first and second terms of Eq. (|2|). After per- 
forming thermodynamic averages one obtains 



In 2 , ^ N , / sinh £ 
— -ln(25+l)+ln(^^^ 

^ 1252 



S 
f O 



(40) 



for ferro- and antiferromagnets, where B = coth^ ~ 1/C 
is the Langevin function. This formula was obtained ear- 
lier 1^ with the help of the Wigner-Kirkwood expansion 
which avoids using effective classical Hamiltonians. It 
strongly resembles the result for the two-spin model, Eq. 
(20) of Ref. d), where the factor 7 is replaced by 6. For 
the energy per spin U = —d\nZ/d{PN) one obtains 



U 
J 



{U^B + 9){B 



52 I 
iB') T {2B 



-eS')]- (41) 



Let us consider now the energy U per spin which fol- 
lows from linear spin-wave theory 



U 



Jo Jo 



2 2S J (27r) 
/3 = l3Jo/S, = iS/Jo)sq 



{271^ exp(/3eq) - 1 



1-A, 



q' 



'1-A2, AF. 



(42) 



Here Jq is the zeroth Fourier component of and 
Aq = Jq/Jfl. The first term in the expression for U is 
the classical ground-state energy, the second term is the 
quantum correction to the former, and the last term is the 
temperature-dependent magnon contribution. Strictly 
speaking, linear spin-wave theory is only applicable for 
higher then two dimensions, but, although, Eq. (^) re- 
mains well-defined in lower dimensions. The Haldane 
gap in the magnon spectrum for integer 5, which is not 
taken into account in Eq. (^), behaves as exp(— 5) and 
becomes negligible for large spins. Whereas Eq is of order 
unity, the parameter /3 is precisely the small parameter 
of the cumulant expansion. If one expands U in pow- 
ers of (3, one obtains a series which is very close to the 
one following from the cumulant expansion in the limit 
^ = J/T ^ 1. In particular, for a ferromagnetic chain in 
the temperature interval J/5' < T < J Eq. (^) yields 



U 
J 



1 

s 



252 



1 

252' 



(43) 



whereas expanding Eq. ( |4^ ) we obtain the same expres- 
sion without the last term. In the antiferromagnetic case 
one obtains similar expressions with coefficients 6 and 3 
instead of 2 and 2 in the denominator. 

In Fig. H we compare the heat capacity C — dU jdT of 
a Heisenberg antiferromagnetic chain with S = "ij^ cal- 
culated in different orders of the cumulant expansion [see 



6 



c 



1 2 -| I I ' I I I ' I ' L 




t= TI{JS^) 



FIG. 3. Heat capacity of the antiferromagnetic Heisen- 
berg 5* = 3/2 chain. Theoretical results to zeroth, first and 
second order of the cumulant expansion, as well as of lin- 
ear spin-wave theory. They are compared with the numerical 
TMRG calculation of Ref. 0. 



physical quantities as the Wigner-Kirkwood expansion, 
whereas the approach of Refs. Q is nonperturbative. 

The effective classical Hamilton function TL obtained 
in this paper can be used to compute quantum correc- 
tions for magnetic systems with large spins and not too 
low temperatures. For one-dimensional models, this was 
done with the help of analytical methods. In two and 
higher dimensions, one can apply the diagram technique 
for classical spins (see, e.g., Ref. Q), which should 
be generalized to non-Heisenberg Hamiltonians, how- 
ever. For models without long-range order, such as two- 
dimensional ferro- and antiferromagnets, it is problem- 
atic to sum up the relevant diagrams for the correspond- 
ing classical case when the temperature is low. Then the 
1/-D expansion, where D is the number of spin compo- 
nents, proves to be an efficient tool for low-dimensional 
classical magnets p^-p^] . Since TL is written in terms of 
various scalar products, it can be easily generalized to 
arbitrary values of D and treated with the help of a 1 /D 
expansion. 



Eq. (^)] and following from the linear spin-wave theory 
[see Eq. (42)] with numerical result of Ref. where 
the transfer-matrix renormalization group (TMRG) was 
used. Although a spin value 3/2 is not large, one no- 
tices that taking into account quantum corrections im- 
proves the behavior of the heat capacity provided the 
temperature is not too low. In order to achieve quantita- 
tive agreement with the numerically accurate result, one 
must take into account higher-order corrections in \/S. 
This requires calculation of the next terms of the effective 
classical Hamiltonian, which can be done with the help 
of the diagram technique developed in Sec. III. Also on 
the low-temperature side the accuracy of the SWT result 
can be improved by taking in Eq. ( [42| ) into account \/S 
corrections to the magnon spectrum. 



VII. DISCUSSION 

In the preceding part of the paper, we have obtained an 
effective classical Hamiltonian H for the large-spin quan- 
tum magnet described by the Heisenberg Hamiltonian 
H, Eq. (|l^). Ti consists of a purely classical part and 
quantum corrections of different orders in 1/5. Quan- 
tum corrections have a non-Heisenberg form and their 
structure becomes more and more complicated with in- 
creasing order of 1/5. In particular, pair interactions in 
H give rise to many-spin interactions in Ti. 

Our effective classical Hamilton function Ti looks dif- 
ferent from the effective classical Hamilton function ob- 
tained in Refs. The latter has the same Heisenberg 
form, and the quantum effects are absorbed in the fac- 
tor renormalizing the exchange interaction. This factor 
satisfies a system of nonlinear equations. It is difhcult 
to make a direct comparison of the two effective classical 
Hamilton functions because of their different structures 
and different ways of derivation: Our approach is per- 
turbative in 1/5 and it leads to the same results for the 



APPENDIX: ELIMINATION OF TRANSVERSE 
COMPONENTS OF THE COHERENT-STATE 
VECTORS 

There are two generic formulas for the elimination of 
transverse components of the coherent-state vectors. The 
first one, 

^ (a • n„)(n„ • b) - (a • b) - (a • n)(n • b), (44) 



where a and b are arbitrary vectors, follows from the def- 
inition of the scalar product (a • b) . The second formula 
is 

(a • nj:){ny • b) - (a • nj^)(n^ • b) = ([n^; x iij^] ■ [a x b]) 
= (n • [a X b]). (45) 

Combining these two formulas one obtains the relation 

4(a • n±)(nzp • b) = (a • b) - (a • n)(n • b) 

±i(n-[axb]), (46) 

which is used in the main text to eliminate n±. Other 
useful relations are 

(n+ -n^) = ^{nl + nD = ^ 

[n+ X n_] = -n, [n± x n] = =F«n±. (47) 

* www.mpipks-dresden.mpg.de/~garanin/; 

garanin@mpipks-dresden.mpg.de 
^ cnIs.lanl.gov/~kladko/; kladko@lanl.gov 
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